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Abstract 

In this article, we take the point of view that the charmed mesons D s o(2317) 
and D s i(2460) with the spin-parity + and 1 + respectively are the conven- 
tional cs states, and calculate the strong coupling constants Gs( for (D*<j)\D s o) 
) and Ga{ for {D s (p\D s i) ) in the framework of the light-cone QCD sum rules 
approach. The strong coupling constants Gs and Ga are related to the basic 
parameter fi in the heavy quark effective Lagrangian, the numerical value is 
larger than the existing estimation. With the assumption of the vector me- 
son dominance of the intermediate 0(1020), we study the radiative decays 
D S Q — > D*7 and D s \ — > Dg'j. 
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1 Introduction 

The observation of the two charmed resonances _D s o(2317) in the D s 7r° invariant 
mass distribution and D s i(2A60) in the D*n° and Z) s 7 mass distributions has trig- 
gered hot debate on their nature, under-structures and whether it is necessary to 
introduce the exotic states [H [2]. They can not be comfortably identified as the 
quark-antiquark bound states in the spectrum of the constituent quark models, 
their masses are significantly lower than the values of the + and 1 + states respec- 
tively from the quark models and lattice simulations [3]. The difficulties to identify 
the _D s0 (2317) and _D s i(2460) states with the conventional cs mesons are rather sim- 
ilar to those appearing in the light scalar mesons below lGeV. The light scalar 
mesons are the subject of an intense and continuous controversy in clarifying the 
hadron spectroscopy [I], the more elusive things are the constituent structures of 
the /o(980) and ao(980) mesons with almost the degenerate masses. The mesons 
D s o(2317) and _D s i(2460) lie just below the DK and D*K threshold respectively, 
which are analogous to the situation that the scalar mesons ao(980) and /o(980) 
lie just below the KK threshold and couple strongly to the nearby channels. The 
mechanism responsible for the low-mass charmed mesons may be the same as the 
light scalar nonet mesons, the / (600), / (980), a (980) and K*(800) [5j [6j El E] . 
There have been a lot of explanations for their nature, for example, the conventional 
cs states [ID] , two- meson molecular states [H], four-quark states [12] . etcH. If we 
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2 The literatures listed here are far from complete, for more literatures, one can consult Ref. [2]. 



take the scalar mesons a (980) and / (980) as four quark states with the constituents 
of scalar diquark-antidiquark sub-structures, the masses of the scalar nonet mesons 
below lGeV can be naturally explained [HIE]- 

There are other possibilities besides the four-quark state explanations, for ex- 
ample, the scalar mesons ao(980), /o(980), D s0 (2317) and the axial- vector meson 
_D s i(2460) may have bare P— wave qq and cs kernels with strong coupling to the 
nearby thresholds respectively, the S— wave virtual intermediate hadronic states (or 
the virtual mesons loops) play a crucial role in the composition of those bound states 
(or resonances due to the masses below or above the thresholds). The hadronic dress- 
ing mechanism (or unitarized quark models) takes the point of view that the mesons 
/o(980), a (980), D s0 (2317) and _D s i(2460) have small qq and cs kernels of the typical 
qq and cs mesons size respectively. The strong couplings to the virtual intermediate 
hadronic states (or the virtual mesons loops) may result in smaller masses than the 
conventional scalar qq and cs mesons in the constituent quark models, enrich the 
pure qq and cs states with other components [131 [T4] . Those mesons may spend part 
(or most part) of their lifetime as virtual KK, DK and D*K states [31 l6l \7\ [13| [T^] . 

The radiative decays can be used to probe the under-structures of the hadrons, 
and they are suitable to understand the nature of the D s0 (2317), -D s i(2460) and 
distinguish among different interpretations JT3J [HH [T71 US EE 120] ■ Different under- 
structures can lead to different decay widths, and the predictions can be compared 
with the experimental measurements. For example, the value of the strong coupling 
constant Qd s qDK with the assumption that the D s0 (2317) being a conventional scalar 
cs state is much larger than (or several times as large as) the corresponding value 
with the assumption of being a tetraquark state [01 [21] ■ The D s0 (2317) can not 
decay to the D s ^j due to the angular momentum and parity conservation, and such 
a final state has not been observed; the decay D s q — > D*^ is allowed and no evidence 
is reported yet of the final state -D s 77 resulting from the decay chain D s0 — > D*^ — > 
-D s 77. The radiative decay widths of the D s0 — > D*^ and D sl — > D s ^ have been 
calculated with the constituent quark model [TT1 [TSJ [TH] , the vector meson dominance 
(VMD) ansatz [16] and the light cone QCD sum rules [13] . etc. 

The amplitudes of the radiative decays D s0 — > D*j and D s \ — ► D s j can be 
written as 



due to the Lorentz covariance. The p^ and k^ are the four momenta of the D*(D s i) 
and 7, respectively; the rj^ and e M are the polarization vectors of the Dl(D s i) and 
7, respectively. The parameters ds and d,A have the dimension of inverse of the 
mass, and get contributions from the photon couplings both to the light quark 
part c s s^^s and to the heavy quark part e c c7 M c of the electromagnetic current, 
here the e s and e c are strange and charm quark charges in units of e. In order to 
determine the amplitudes of the D s0 — > D*^ and D s \ — > Dgj, we follow the VMD 
ansatz [161 122]- In the heavy quark limit, the matrix elements (D*(v', rf)\cq s q{v)) 



(D* s (p, 77)7(^,6)1^0) 
(D a7 {k,e)\D 8l (p,Ti)) 



ed s {e* ■ r)*p 
iedA {e* • rjp 



k — e* ■ prf ■ k} , 
k — e* • prj ■ k} 



(1) 
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and (D s (v') | c^^c\D sl {v,7])) vanish for v ■ v' = (M^ s0 + Mf,*) /2M Ds0 M D * w 1 and 
v ■ v' — {Ml) si + Mj) a ) /2M Dsl M Ds « 1, here the v^, v' are the four-velocities of the 
heavy mesons. We consider only the contribution of the intermediate 0(1020), 

(^(p^MMIAk)) = WM<Kl> 0I^)^<7(*, e)|0(g,O) 



2 

(.P;(p, ry)0(fc, 01 Ao) 7 2 _ m 2 U m ^Q s {-i)e* ■ £ 



= G s (e* - V *p-k-e* -pn* -k)^eQ s , (2) 

= ( D s<l>{k,£,)\D sl (p,r)))— ^f^eQ^-^e* •£ 

= iG A (e* ■ r)p- k - e* ■ prj ■ k)—eQ s , (3) 

where we have use the VMD Lagrangian 

Ml 

C = -eQ—A a <p a = -eQUM^r , (4) 
and the definitions of the strong coupling constants 

(D'sMtMlDa) = G s (C-v*p-k-e-Pn*-k), 

(D s <p(k,0\D sl (p, V )) = iG A (C-VP-k-C-PV-k), 
(O|s(O) 7m s(O)|0(A:,O> = Um&- ( 5 ) 

The G5 and G A are the strong coupling constants, the is the weak decay constant 
of the vector meson <p and the £ M is the polarization vector. The strong coupling 
constants Gs and G A can be related to the effective coupling constant fi in the 
heavy quark effective Lagrangian. The Lagrange density is set up by the hidden 
gauge symmetry approach with the light vector mesons collected in a 3 x 3 matrix 
% [22] , 

C' = ijlTr{S a H b a at3 V^} , (6) 

where the effective fields H a and S a stand for the doublets with J p = (0~, 1~) and 
(0 + , 1 + ) respectively, 



Sa = ^ [Jftv* - ^a] , (7) 
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where the v M is the four- velocity of the heavy meson and the a is a light quark flavor 
index. The P 0a , P a , P* and P[ a are the scalar, pseudoscalar, vector and axial-vector 
mesons respectively The H a = ^H^ , S a = 7 ^7°, V a p = d a V - d p V a + [V a , V ] 
and V a = i^V a , the gy is fixed to be gy = 5.8 by the KSRF rule [23]. Finally we 
obtain the relation between the G$ {Ga) and fx 

G s = 2V2figv, 

G A = 2V2ftg v . (8) 

The parameter fx is a basic parameter in the heavy quark effective Lagrangian, the 
precise value can lead to more deep understanding of the relevant dynamics. It 
is interesting to calculate its value with the light-cone QCD sum rules. The fx is 
estimated to fx = — 0.13 ±0.05 GeV^ 1 from the analysis of the D — > K* semileptonic 
transitions induced by the axial weak current [22J. 

In this article, we take the point of view that charmed mesons D s0 (2317) and 
_D sl (2460) are the conventional cs states, calculate the values of the strong coupling 
constants Gs and Ga ( and the corresponding fx ) in the framework of the light- 
cone QCD sum rules approach, and study the radiative decay widths of the D s0 — > 
D*-f and D sl -> D s ^.. The light-cone QCD sum rules approach carries out the 
operator product expansion near the light-cone x 2 m instead of the short distance 
x ~ while the non-perturbative matrix elements are parameterized by the light- 
cone distribution amplitudes which classified according to their twists instead of the 
vacuum condensates [2U [25]. The non-perturbative parameters in the light-cone 
distribution amplitudes are calculated by the conventional QCD sum rules and the 
values are universal |26j . 

The article is arranged as: in Section 2, we derive the strong coupling constants 
Gs and Ga within the framework of the light-cone QCD sum rules approach; in 
Section 3, the numerical result and discussion; and in Section 4, conclusion. 



3 The same approach can be used to study the radiative decay widths of the D s i — > D*j, 
D s i — > D s0 -f, and explore the structures of the mesons D s i(2460), D s0 (2317), that may be our 
next work. 
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2 Strong coupling constants Gs and Ga with light- 
cone QCD sum rules 

In the following, we write down the two-point correlation functions U^(p,q) and 

n>,g), 

KM = ijd'xe-^iOlTlJ^O^ix)}^)), (9) 

TLfaq) = iJfae^MTiJtWMz)}^)), (10) 

J (x) = s(x)c(x) , 
J 5 (x) = c(x)ij 5 s(x) , 



J»{x) = c{x) ltl s{x) 



tf(x) = s{x) W c{x), (11) 

where the currents J s (x), J$(x), JY(x) and Jfi(x) interpolate the mesons D s0 (2317), 
D s , D* and £) s i(2460), respectively, the external state <p has the four momentum p^ 
with p 2 = m|. The correlation functions IH (p, q) and n^(p, g) can be decomposed 
as 

nj(p,g) = n 5 {e M g • p - p M e • g} + ••• , 

nj(p,g) = in A {e M (g + p).p-p M e.(g + p)} + ... (12) 

due to the Lorentz covariance, here the e M is the polarization vector of the <f) meson. 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach [26J, we can insert a complete series of intermediate states with the 
same quantum numbers as the current operators J s (x), J§(x), JY(x) and Jf}(x) into 
the correlation functions EP (p, q) and n^(p) to obtain the hadronic representation. 
After isolating the ground state contributions from the pole terms of the mesons 
D s0 (2317), D s , D* s and D s i(2460), we get the following results, 

" P,(1 {M 2 Ds0 -{q + P?}{Ml t -q 2 ) 

= f^M Da0 M Dt G s { .^..^ (13) 

{Ml sQ -{q+pY}{Ml t -q>) 

uA( , (0\J*(0)\Dsi(q + p))(Dsi\Ds<f>)(Ds(q)\J 5 (0)\0) 

AP,(1) {Ml sl -{q + p?}{Ml s -q 2 } 

= l l I ^ Tm^~T*T^J a#2 txMl + P) -P-Ppe- (q + p)} 

(m c + m s ) [M z Dsi - (q + p) 2 } [M 2 Ds -q 2 \ 

+ -.., (14) 
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where the following definitions have been used, 



(0\J- 



(0)\DM) 



fD s0 M Ds0 , 

f D M 2 Ds 
m c + m s ' 



<0|J B + (0)|£>.(p)> 
(0\J*(0)\D sl (p, V )) 



fD sl M Dsl 7]^ 



where the /d s0 , /d s , /d» and /x> al are the weak decay constants of the D s0 (2317), D s , 
D* and _D s i(2460), respectively. In Eqs.( 13-14), we have not shown the contributions 
from the high resonances and continuum states explicitly as they are suppressed due 
to the double Borel transformation. 

In the following, we briefly outline the operator product expansion for the corre- 
lation functions U.^(p, q) and n^(p, q) in perturbative QCD theory. The calculations 
are performed at the large space-like momentum regions (g + p) 2 <C and q 2 -C 0, 
which correspond to the small light-cone distance x 2 ~ required by the validity 
of the operator product expansion approach. We write down the propagator of a 
massive quark in the external gluon field in the Fock-Schwinger gauge firstly [27] . 



here we have neglected the contributions from the gluons G^ v . The contributions 
proportional to the G^ u can give rise to three-particle (and four-particle) meson 
distribution amplitudes with a gluon (or quark- ant iquark pair) in addition to the 
two valence quarks, their corrections are usually not expected to play any signif- 
icant rolefl Substituting the above c quark propagator and the corresponding 
meson light-cone distribution amplitudes into the correlation functions IT^Qo, q) and 
rr4(p, q) in Eqs.(9-10) and completing the integrals over the variables x and k, finally 

4 For examples, in the decay B — > \coK, the factorizable contribution is zero and the non- 
factorizable contributions from the soft hadronic matrix elements are too small to accommodate 
the experimental data |28j ; the net contributions from the three- valence particle light-cone distri- 
bution amplitudes to the strong coupling constant (Jd s iD*k are rather small, about 20% [7|. The 
contributions of the three-particle (quark-antiquark-gluon) distribution amplitudes of the mesons 
are always of minor importance comparing with the two-particle (quark-antiquark) distribution 
amplitudes in the light-cone QCD sum rules. In our previous work, we study the four form-factors 
fi{Q 2 )i f2{Q 2 ), 9i{Q 2 ) and g2(Q 2 ) of the S — * n in the framework of the light-cone QCD sum 
rules approach up to twist-6 three-quark light-cone distribution amplitudes and obtain satisfactory 
results [29]. In the light-cone QCD sum rules, we can neglect the contributions from the valence 
gluons and make relatively rough estimations. 




(15) 
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we obtain the results 
■l 



n 



s 



+ 



7 



4 

(T2 



AA 



fl 



duA i (u) 



1 



+ 



2m s 
2ml 



m„ 



-fimi / duu / dtC ± (t) 



AA 2 AA 3 
1 



+ 2/>, 



y 2 ' ' du 







+ 



2m 



2 i 



+2f (f) m (j) m c / dwu 



AA 2 AA 3 
A A 3 ' 



- 2f <j) m 4> m c 



dr 
i 



duu 



o 



AA 2 

(«) 

AA 2 



(16) 



n 



where 



o 



T^2 />! 



4 

fT 2 



(IwA i (u) 



f<t> f<t> 

1 



+ 



m<t> 
2ml 







/^mj / dwx / dtC ± (t) 



AA 2 AA 3 
1 



AA 2 
i 



+ 2/> 2 I du dr 

'o </o 
• i 







gig) 
AA 2 



+ 



2m 2 



+2f ( f ) m ( j > m c I duu I dr dt 
Jo Jo Jo 



AA 2 AA 3 



2f <t> m !t) m c 



duu 



(v) 

9± 
AA 2 



AA 3 



(17) 



AA 



m c — (q + up) . 



In calculation, the two-particle meson light-cone distribution amplitudes have 
been used [30], the explicitly expressions are given in the appendix. The parameters 
in the light-cone distribution amplitudes are scale dependent and can be estimated 
with the QCD sum rules approach [30] • In this article, the energy scale fi is chosen 
to be fi — lGeV. 

Now we perform the double Borel transformation with respect to the variables 
Qi = ~(P + <?) 2 an d Q\ = —<f for the correlation functions lis and ILa in Eqs.(13- 
14), and obtain the analytical expressions of the invariant functions in the hadronic 
representation, 



B M 2B A/l 2li 



M 1 2ii 5 



G s f Dt f Da0 M Dl M DsQ 

M'rMj 

G A fpjD sl M 2 D M D sl 
MlMl{m c + m s 



exp 



exp 



Ml 


Ml, 


Ml 


M| 


r mi 


MIA 


Ml 


M| J 



+ ■ 



+ 



(19) 



here we have not shown the contributions from the high resonances and continuum 
states explicitly for simplicity. In order to match the duality regions below the 
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thresholds s and s' Q for the interpolating currents J (x) (or j£(x)) and Jj(x)(or 
Js(x)) respectively, we can express the correlation functions II (denote the Us and 
Ha) at the level of quark-gluon degrees of freedom into the following form, 

ds [ ds'- ^ 



n - J"* J "jrtoWw^Ty (20) 

where the p(s, s') are spectral densities, then perform the double Borel transfor- 
mation with respect to the variables Q\ and Q\ directly. However, the analytical 
expressions of the spectral densities p(s, s') are hard to obtain, we have to resort 
to some approximations. As the contributions from the higher twist terms are 
suppressed by more powers of m 2_(^ +up y , the net contributions of the twist-3 and 
twist-4 terms are of minor importance, less than 20%, the continuum subtractions 
will not affect the results remarkably. The dominating contribution comes from 
the two-particle twist- 2 term involving the 4>±(u). We preform the same trick as 
Refs. [27| I3T] and expand the amplitude <j)±(u) in terms of polynomials of 1 — it, 

± (n) = ^6 fc (l-u) fc = ^6 fc (^^ , (21) 

fc=0 k=0 ^ S ^ ' 

then introduce the variable s' and the spectral density is obtained. 

After straightforward calculations, we obtain the final expressions of the dou- 
ble Borel transformed correlation functions n(Mf,M|) at the level of quark-gluon 
degrees of freedom. The masses of the charmed mesons are Mr> sl = 2A6GeV, 



M Ds0 



2.317GeV, M Da = 1.97GeV and M D * = I.WIGeV 



Mr 



' M D *+M D . 



0.48, 



0.45, there exists an overlapping working window for the two Borel 



M Ds +M Dal 

parameters Ml and Mf , it's convenient to take the value Mf 
the threshold parameter so and make the simple replacement, 



Mf . We introduce 



m c+"0( 1 - u 0) m </ ) 



m c+ u 0( 1 - u 0) m 1 



s 
A/ 2 



to subtract the contributions from the high resonances and continuum states 
finally we obtain the sum rules for the strong coupling constants G$ and Ga, 



G s f Ds0 f Da M Ds0 M D *exp 



fTM 2 <j> ± (u ) <|exp 



exp 



BB 



fl 



BB 

iM 2 

„ 2m 



M 2 n 


Ml, 


Mf 


Mf 




r s 


— exp 


I M 2 



f$mlA ± (u ) 



1 + 



M 2 



+ 



pu Q pT pu 

/ dr dtB ± (t) - /Jmj^ / dtC ± (t) 
Jo Jo Jo 

o M 2 j 



m; 
1 H - 

M 2 



2f 4> m (t) 'm c uog]_ > (uo) + f^m^rricUo 



dr 



(22) 
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GAfD al fD s M Dgl 



M 2 Ds 



m c + m. 



exp 



M 2 Dsl Ml a 



M{ 



Ml 



-f'M z Muo)\exp 
BB 



BB 



exp 



s 



M 2 



+ exp 



M 2 

Uo 



rp 2lTl S 

m <t> 



f$m 2 A ± (u ) 



1 + 



+2^mJ I dr J dtB ± (t) - fimfa \ dtC ± (t) 

cm 

W 



i + 



rn 



2 i 



M 2 



dr 



where 



M 2 



(23) 



BB 

u 

M 2 



m 2 c +u (l -u )ml 
M 2 



Mf + Mf ' 

M 2 M 2 
M 2 + Mf ' 



(24) 



3 Numerical result and discussion 

The parameters are taken as m s = (140 ± 10)MeV, m c = (1.25 ± 0.10)GeV, = 
(0.215 ± 0.005)GeV, f£ = (0.186 ± 0.009)GW, a£ = 0.2 ± 0.2, a| = 0.2 ± 0.1, 
q 3 = 0.032 ± 0.010, ^ = 0.15 ± 0.10, = 0.10 ± 0.05, <f = -0.10 ± 0.05, = 
-2.1 ± 1.0, ujV = 3.8 ± 1.8, to? = 7.0 ± 7.0, m = 1.02GeV, M Ds0 = 2.317GeV, 
M Dsl = 2A6GeV, M D * = 2.112GeV, M Ds = 1.97GeV, f D . = (0.26 ± 0.02)GeV, 
f Ds = (0.26±0.02)GeV, f Ds0 = (0.225 ± 0.025)GW and f D J x = (0.225 ± 0.025)GW 
[IH [32]. The duality thresholds s in Eqs. (22-23) are taken as s = (7.0 - 8.0)GeV 2 
(^/so = (2.6 — 2.8)GeV) to avoid possible contaminations from the high resonances 
and continuum states, in this region, the numerical results are not sensitive to 
the threshold parameters Sq. The Borel parameters are chosen as Mf = Mf and 
M 2 = (3.5 — 8)GeV 2 , in those regions, the values of the strong coupling constants 
Gs and Ga are rather stable from the sum rules in Eqs. (22-23) with the simple 
subtraction. 

The uncertainties of the eleven parameters f^, , a|, uj, , ui^, <j 3 , q, <^J, 
and m s can not result in large uncertainties for the numerical values. The main 
uncertainties come from the six parameters /d s0 , /d s1 , /dj, /d s , ^ and m c , the 
variations of those parameters can lead to relatively large changes for the numerical 
values, while the dominating uncertainty comes from the a^, which are shown in 
the Fig.l, refining the parameter is of great importance. Taking into account 
all the uncertainties, finally we obtain the numerical results of the strong coupling 
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Figure 1: The Gs and Ga with the Borel parameter M 2 for = 0.2 (A), = 0.4 
(B) and 4 = 0.0 (C). 

constants Gs and G^, which are shown in the Fig.2, 

1 



G s 
G A 



-(4.4 ±2.0) 
-(3.9 ±1.6) 



GeV 
1 

GeV 



-(0.27 ±0.12) 
-(0.24 ±0.10) 



GeV 
1 

GeV 



(25) 



(26) 



The corresponding values of the parameter fx are larger than the existing esti- 
mation fx = —0.13 ± 0.05 GeV" 1 from the analysis of the D — > K* semileptonic 
transitions induced by the axial weak current [22J . From the numerical values of the 
strong coupling constants Gs and Ga, we can obtain the decay widths, 



T(D s0 -> D* sl ) 
T(D sl -> D sl ) 



{1.3 - 9.9)KeV , 
(5.5-31.2).KeV. 



(27) 



The comparison with the results from other approaches is presented in the Table. 1. 



4 Conclusion 

In this article, we take the point of view that the charmed mesons _D s o(2317) and 
-D s i (2460) with the spin-parity + and 1 + respectively are the conventional cs states 
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Figure 2: The Gs and with the Borel parameter M 2 , the A, B and C stand for 
the central value, down limit and up limit, respectively. 



KeV 




PS] 


[16J 


[17j 


[18j 


[19J 


[20j 


This work 


T(D s0 - £> s * 7 ) 


4 


-6 


1 


1.74 


1.9 


1.1 


35 


1.3-9.9 


r(D sl -> D sl ) 


19 


- 29 




5.08 


6.2 


0.6-2.9 




5.5-31.2 



Table 1: The decay widths of the D s0 — > D*j, D si — > D s ^. 



and calculate the strong coupling constants Gs and Ga in the framework of the 
light-cone QCD sum rules approach. The strong coupling constants Gs and Ga 
are related to the basic parameter fx in the heavy quark effective Lagrangian, the 
numerical value of the fx is larger than the existing estimation. With the assumption 
of the vector meson dominance of the intermediate 0, we study the radiative decays 
D s0 — > Z)*7 and D s i — > D s j, the numerical values of the decay widths are compatible 
with the existing estimations, further experimental data can conform or reject the 
assumption of the two-quark substructure. Just like the scalar mesons /o(980) and 
ao(980), the scalar meson D s q(2317) and the axial-vector meson Z) s i(2460) may have 
small cs kernels of typical cs meson size. The strong couplings to virtual intermediate 
hadronic states (or the virtual mesons loops) can result in smaller masses than the 
conventional + and 1 + mesons in the constituent quark models, enrich the pure cs 
states with other components. The -D s o(2317) and D sl (2A60) may spend part (or 
most part) of their lifetimes as virtual DK and D*K states, respectively. 
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Appendix 

The light-cone distribution amplitudes of the meson are defined by 

(0\s(0)^„s(x)\(f)(p)) = p^m^ 



(O|s(O)s(x)|0(p)) 



+ 



1 



e ■ x 



p ■ x 
1 e • x 



Urn* / due- lupx gl>{u) 



o 



— -x 



r^.^l *«-*-C(«), 

j, 2,TYl s 



m\e-x I due- lup - x h\°\u), 



o 



-l 

fT 2 



2 2 

mix 



e • x 



- p^]/^ "l^^T^a 



due- lup - x B ± (u) 



o 



-tup-x 



C ± (u). (28) 
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The light-cone distribution amplitudes are parameterized as 
^(u,^ = 6u(l-u){l + a 2 ^(5£ 2 -l)j , 

6u(l-u){l + 4^(5^-1)} , 
^(l + a + {^l + 5. 3 }(3e 2 -l) 

+ {^4 + ^ft(3w s v " ^)} (3 - m 2 + 35^ 4 ) , 

f 2 H 40 20 1 i 

g 3 (u,(i) = l + |-l-^ + y?3-yQ|C 2 2 (0 

h 3 (u,fx) = l + |-l + ^a^-10( ?4 T + ?f)}c|(O 

/ij'W) = 6 M (l- M )|l + (ia| + ^ 3 ^3 T )(5e 2 -l)} , 
JifW) = 3e 2 + ^ 2 (5e 2 -3) + ^ 3 T (3-30e 2 + 35e 4 ), 

A ± (u,fi) = 30u 2 (l-u) 2 || + ^4 + t% ~\%\ , 
C(u,fi) = g 3 (u,fj) + 4>\\(u,fi) - 2g ( f\u,fi) , 

B±(u,fi) = h ( l\u,fi) - -<t>±{u,fj,) - -h 3 (u,fi) , 

C±(u,fx) = h 3 (u, //) - 4>±(u, //) , (29) 

i j_ 

where the £ = 2it — 1, and the C 2 2 and C| are Gegenbauer polynomials [30J. 
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